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A 20-kiloton liquid scintillator detector is designed in the Jiangmen Underground Neutrino Observatory 
(JUNO) for multiple physics purposes, including the determination of the neutrino mass ordering through reac- 
tor neutrinos, as well as measuring supernova neutrinos, solar neutrinos, and atmosphere neutrinos to explore 
different physics topics. Efficient reconstruction algorithms are needed to achieve these physics goals in a wide 
energy range from MeV to GeV. In this paper, we present a novel method for reconstructing the energy of 
events using hit information from 3-inch photomultiplier tubes (PMTs) and the OCCUPANCY method. Our 
algorithm exhibits good performance in accurate energy reconstruction, validated with electron Monte Carlo 
samples spanning kinetic energies from 10 MeV to 1 GeV. 
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I. INTRODUCTION 


Liquid scintillator (LS) detectors are extensively used in 
nuclear physics and particle physics. Over the past few 
decades, LS detectors have played crucial roles in achieving 
remarkable scientific results in neutrino experiments [1-5]. 
The central detector (CD) of Jiangmen Underground Neutrino 
Observatory (JUNO) will be the world’s largest liquid scin- 
tillator detector, aiming to probe multiple physics goals in- 
cluding determining the mass ordering of neutrinos and accu- 
rately measuring neutrino oscillation parameters through re- 
actor antineutrinos, as well as observing supernova neutrinos, 
solar neutrinos, atmospheric neutrinos, etc [6]. JUNO has 
developed highly transparent LS and highly efficient PMTs 
with 78% photo coverage. These optimizations are necessary 
for JUNO to achieve the key requirements for determining 
the neutrino mass ordering, including an unprecedented en- 
ergy resolution of 3%/./£ (MeV) and better than 1% energy 
scale uncertainty. Additionally, JUNO has developed a dual 
calorimetry technique [7, 8] that can not only calibrate the 
non-linearity of the charge response of 20-inch PMTs, but 
also to enable the detector to operate over a larger dynamic 
energy range. 

Meanwhile, efficient algorithms are necessary for the re- 
construction of individual event energies in JUNO. For the 
energy region of the inverse beta decay (IBD) of reactor neu- 
trinos (Eyis < 10 MeV), JUNO has developed many robust 
reconstruction algorithms based on traditional methods [9- 
11] or the machine learning technology [12, 13]. Events with 
visible energies larger than several hundreds of MeV in the 
detector are primarily from atmospheric neutrino interactions, 
and their energy and direction reconstructions have been well 
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studied using both the probabilistic unfolding method [14, 15] 
and the machine learning technology [16, 17] to assist in de- 
termining the neutrino mass ordering [18]. However, energy 
reconstruction in the mid-energy region (10 MeV < Eķyis < 
several hundreds of MeV), which includes events from dif- 
fuse supernova neutrino background (DSNB) [19], Michel 
electrons, low-energy atmospehric neutrinos, and possible 
proton decays [20], is rarely discussed in previous studies. 
The basic idea of current energy reconstruction algorithms 
in the JUNO central detector are based on either a data-driven 
maximum likelihood method, or a machine learning strategy 
that utilizes the information from the detector hit pattern. The 
data-driven maximum likelihood method has the advantage 
of better modeling the response of the real detector. It pri- 
marily utilizes the detector response to radioactive sources as 
calibration templates, with the energy of radioactive sources 
being in the MeV range [23]. However, it is found that the ac- 
curate energy reconstruction below 10 MeV can be affected 
by the spatial scale of energy deposition [11]. Therefore, it is 
necessary to investigate the feasibility and performance when 
applying this method to cases where events cannot be treated 
as point-like. We introduce the second moment S (Eq. 1) to 
describe the shape of both point-like and cluster-like events. 
This physical quantity is commonly used in accelerator exper- 
iments to describe the shape of clusters in energy calorimetry. 


Gy AE X [re (£a, Yas Za) = ? (a, Y, z)? 
Laci ARa 

where Na g is the number of secondary energy deposition 

for the event. AE, and 72(@a, Yo; Ža) are the energy depo- 

sition and position in the at” secondary energy deposition, 

respectively. P(x, y, z) is the energy-deposit center for the 


event, which is the weighted average of secondary energy de- 
position and can be calculated as follows. 


S= 


(1) 


TONAR (£a Yas Za) X AEa 


Na 
piers AE 

Figure | compares the distributions of the second moment 
for electrons with different kinetic energies in the JUNO’s 


(x,y,z) = (2) 


67 LS. As the kinetic energy of the electron increases, the dis- 
es tribution of the second moment becomes more diffuse, which 
indicates that it corresponds to a larger cluster of energy depo- 
sition. For comparison, Fig. | also shows the second moment 
distribution of muons with different energies deposited in the 
LS. It can be found that the second moment distribution of 
the cluster-like events and the track-like events with the same 
energy deposition is very different, which is related to their 
shape differences in energy deposition. Moreover, for energy 
reconstruction of high-energy events, the potential deviation 
of PMT’s reconstructed charge (charge non-linearity) is also 
an important issue to be considered. 
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Fig. 1. Distributions of the second moment for electrons with differ- 
ent kinetic energies 


79 The machine learning approach demonstrates good perfor- 


sæ mance in event reconstruction as it can effectively utilize the 
s1 detector response information. However, it is important to 
s2 note that this method relies on reliable data samples or sim- 
s3 ulation. The former requires reliable selection and data ac- 
s4 cumulation, while the latter usually needs adjustments ac- 
ss cording to real data, especially in the early stages of detector 
ss Operation. And similar to the data-driven maximum likeli- 
s7 hood method, since the PMT hit pattern is the basic input of 
ss Machine learning, investigating the influence from potential 
s charge non-linearity of PMTs based on real data is necessary. 
90 In this paper, we investigate the energy reconstruction of 
ə events spanning from MeV to GeV in the JUNO CD using 
2 the data-driven maximum likelihood method [11]. To re- 
ə duce the dependence on the accuracy and non-linearity of 
ə PMT charge reconstruction, we only utilize information of 
ə PMT firing states (fired or unfired), known as the OCCU- 
os PANCY method. This paper focuses on the reconstruction 
97 Of cluster-like events since the calibration templates in our 
ə Study are constructed using point-like (or cluster-like) cali- 
ə bration sources. Thanks to the study of event identification in 
10 JUNO [14], we can effectively distinguish cluster-like events 
101 from track-like events. As for the reconstruction of track-like 
1o2 events (mainly muon-like events that have long tracks in the 
103 detector), it is also an important topic and has been carried 


1⁄4 Out in [14—-17, 24-27]; however, it is not the subject of this 137 


107 Offline software [28—30] is used for validation. 

The details of our study are presented as follows: First, 
1o9 We introduce the JUNO detector and 3-inch PMT system 
110 (Sec. II). Then, we present the methodology of our recon- 
111 struction (Sec. III), including the construction of calibration 
112 Map and the construction of maximum likelihood function. 
113 In Sec. IV, the reconstruction performance will be shown and 
114 compared. Finally, a summary will be provided in Sec. V. 


108 


Il. JUNO DETECTOR AND 3-INCH PMT SYSTEM 


ne As shown in Fig. 2, JUNO mainly consists of three sub- 
117 detectors: CD, the water cherenkov detector, and the top 
11g tracker detector [6, 31]. CD contains 20 kilotons of LS and 
119 a 12 cm thick acrylic spherical container with a diameter of 
120 35.4 m. The main component of the LS is linear alkyl ben- 
121 zene(LAB), with PPO (2.5-diphenyloxazole) as fluor and bis- 
122 MSB as wavelength shifter. A total of 17612 20-inch PMTs 
123 (LPMTs) and 25600 3-inch PMTs (SPMTs) will be installed 
124 On the exterior of the container as photosensors to collect pho- 
125 ton signals. As a result, more than 1345 photoelectrons (PEs) 
126 Will be observed by CD for a 1 MeV electron that fully de- 
127 posits its kinetic energy in the LS. The SPMTs will work al- 
123 most exclusively in the single photoelectron (spe) mode for 
129 reactor antineutrino detection (Eyis < 10 MeV). Therefore, 
130 SPMTs can serve as a linear reference for LPMTs and be used 
131 to calibrate the charge non-linearity of LPMTs. This feature 
132 is helpful in constraining some of the systematic uncertainties 
133 in the LPMT energy reconstruction and improving the energy 
134 resolution. Moreover, SPMTs have the potential to detect su- 
135 pernova neutrinos and measure the solar parameters (02 and 
130 Am3,) independently [32]. 
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Fig. 2. A schematic view of the JUNO detector [6] 


For the detection of events with energies greater than tens 


105 article. To investigate the performance of our reconstruction, 138 of MeV or even GeV, most LPMTs will receive tens or even 
1 Monte Carlo (MC) simulation data generated by the JUNO 10 hundreds of PEs. For example, Fig. 3(a) shows the distribu- 


14 


So 


14 


142 
143 
144 
14 


a 


14 


a 


14 


à 


14 


© 


14 


© 


151 


So 


15 


152 


15: 


o 


154 


15: 


a 


15 


D 


15 


cil 


15 


© 


15: 


© 


160 


161 


162 


16; 


wo 


164 


165 


16 


D 


16 


=] 


16 


© 


16 


© 


171 


=] 


17 


172 


Detector volume 
all 


— R<10m k 
—— 10m<=R<14m 


— 14m<-R<15m 
—— 15m<-R<16m 
—— Retém 


Detector volume —— 500MeV 
—— 1000MeV 


—— 2000MeV 


oS 
& 
T 


— 14m<-R<15m 
—— 15m<-R<16m 
—— R<i6m 


5 
N 
"Į 


’ 
re 
’ 
’ 
PE diy, 
tii 
tess, 
O4 


pad 
D 


V4 @ 


firedSPMT/allSPMT 


= 

a 
TTT 

p 


Yvy 
Yvy 
vy 
Yvy 


0.46 


fi fi L L fi i I fi L 
O 100 200 300 400 500 600 700 800 900 1000 0 5 15 20 


nPEbyPMT [p.e.] 


(a) nPE received by LPMTs for 500 MeV 
electron 


10 


nPEbyPMT [p.e.] 


(b) nPE received by SPMTs for 500 MeV 
electron 


ni L E 
fi fi fi L fi fi fi fi 
25 30 0.8) 500 1000 1500 2000 , 2500 3000 3500 4000 
R° [m 


(c) The proportion of fired SPMTs 


Fig. 3. The distribution of nPE received by LPMTs (a) and SPMTs (b) for 500 MeV electron deposited its kinetic energy in the LS. (c) The 
proportion of fired SPMTs for electrons deposited their kinetic energies at different locations 


tion of the number of PEs (nPE) received by LPMTs for a 
500 MeV electron that deposits its kinetic energy in the LS. 
Obviously, all LPMTs are fired in this case. If the energy de- 
position occurs at the edge of the detector, the nearby LPMTs 
will receive even more PEs. The linearity of LPMT charge re- 
construction over a large charge dynamic range is a challenge 
and needs to be calibrated and validated in the future. Based 
on the experience of Daya Bay, an independent measurement 
system will be an effective solution [33, 34]. For compar- 
ison, in the same case of a 500 MeV electron deposited its 
kinetic energy in the LS, there is about 45%-60% of SPMTs 
that are not triggered (Fig. 3(c)), and most fired SPMTs only 
receive less than 5 PEs (Fig. 3(b)) due to the fact that their 
photocathode areas are about 40 times smaller. Therefore, we 
develop an energy reconstruction algorithm using only the in- 
formation from SPMTs. In addition, according to the study 
in [35, 36], the readout electronics of JUNO SPMTs may also 
exhibit non-linearity when receiving multiple hits. To mini- 
mize this effect, we will only use its firing information (fired 
or unfired). More details are introduced in Sec. III 


Il. METHOD OF ENERGY RECONSTRUCTION 


A. The probabilities of SPMT’s firing states 


For each SPMT, the number of the detected PEs obeys a 
Poisson distribution: 


eT Hi x pi 


ial (3) 


Poisson(k;|u;) = 


where k; is the nPE detected by the it? SPMT, and li is the 
mean value of the Poisson distribution. As mentioned above, 
in order to reduce the dependence on charge reconstruction 
accuracy, the OCCUPANCY method is applied which uses 
only the information from SMPT’s firing states (fired or un- 
fired). Therefore, only two states of k; need to be considered: 
ki = 0 (unfired) and k; > 0 (fired). The probabilities of these 
two states can be described by Eq. 4 and Eq. 5 [21, 22]. 


Punfirea = Poisson(k; = O|u;) =e” (4) 


173 


Pårea = Poisson(k; > O| ui) 


=1—- Panfired =le“ 


(5) 


174 


In real detection, k; will be smeared by fluctuation in pho- 
toelectron detection, while u; will be distorted due to the ad- 
ditional contribution from PMT dark count. In addition, to 
a avoid false triggering due to electronic noise, we need to ap- 
ply thresholds to each SPMT, which also affect u;. In gen- 
eral, 0.3 PEs is chosen as the typical threshold since it can 
1 effectively handle electronic noise. Taking these effects into 
2 account, the probabilities of unfired state and fired state are: 
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Panfired = P(qi < qitreshal | true) 
= Poisson(ki = Olu) + PihreLoss (pi) 
Parea = P(qi z qld true) 


=1- Punfired 


(6) 


183 


where q; is the reconstructed charge of the it” SPMT, pt” 


184 i 


185 (ute = pPPY + 112”) is the mean value of the Poisson distri- 
186 bution, which consists of two components: 

(1) p? aY caused by the visible energy of physics events; 

188. (2) u introduced by the dark count (DR;) of the ith 
189 SPMT, and it can be calculated by w2" = DR; x t ina time 
190 Window of t; 

191 Pthretoss(ui"*) is the probability of qi 
192 (0.3 PEs in this study) in the case of k; > 0, calculated as 
193 following: 


187 


threshold 
< Gj 


PihreLoss oe) 
i threshold 
= 5 [Poisson(k;| u°) x | Gaus(gi, Cg; ) dq], 
ki=1 0 
with gi = gain X ki, Og = gi x o 


(7) 
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where n indicates the case of multiple PEs, a corre- 
sponds to the ratio between the real SPMT gain and the nor- 


mal SPMT gain (3 x 10°) in JUNO, o;?° denotes SPE resolu- 


tion of the it? SPMT. In real detection, gran, o;P° and DR; 
can be obtained from PMT calibration. 


B. Construction of the calibration map 


JUNO designed a comprehensive calibration system [6] to 
understand the detector response, deploying multiple radioac- 
tive sources in various locations inside/outside of the CD, 
including Auto Calibration Unit (ACU), Cable Loop Sys- 
tem (CLS), Guide Tube Calibration System (GTCS), and Re- 
motely Operated Vehicle (ROV). Figure 4 shows the individ- 
ual calibration systems in the CD and their scanning regions. 
For example, ACU system scans the detector response along 
the central axis with multiple calibration sources, and the CLS 
system can scan in a 2-dimensional plane (X-Z plane) with 
multiple calibration sources using the central cable and side 
cable. The strategy of the JUNO calibration system has been 
developed and optimized based on Monte Carlo simulation 
results [23]. 


source 


we? 


> Water 
Connecting bars 


Fig. 4. The individual calibration systems in the CD and their scan- 
ning regions [10] 


In energy reconstruction, pry directly corresponds to the 
visible energy of an event in the detector, which is the basis 
of the energy reconstruction in our method. Assuming that a 
calibration source is loaded at the location P(r, 0,ġ = 0) in 
the central detector, the mean value of visible-energy-induced 
PEs for the i*” SPMT is ??Y-*°""°°, which corresponds to the 
visible energy (denoted as £*°""°°) of the calibration source. 
Then for an event depositing its energy at the same location, 


the relationship between event’s visible energy Eyis and pry 
for the it” SPMT can be described as follows: 

phy _ vis phy_source 

i -~ ppsource i (8) 


It should be noted that :?"Y is not only related to the visible 


227 energy and the position of the event, but also to the relative 
position (pyr) of the event and the it? SPMT. This rela- 
tionship can be determined using the calibration data, which 
means constructing a calibration map. In Sec. III C, the max- 
imum likelihood method will be adopted to reconstruct the 
visible energy by estimating phy via the calibration map and 
invoking the firing states of 25600 SPMTs from data. Next, 
the construction of the calibration map will be introduced. 

According to the calibration strategy in JUNO, this study 
uses the °*Ge source (E5°% = 1.022 MeV) to calibrate the 
X-Z plane with assistance from both ACU and CLS systems 
across 227 positions. Using the JUNO offline software [28— 
30], 10000 ĉ8Ge events are generated respectively for each 
calibration location on the the X-Z calibration plane. Real- 
istic detector geometry is employed for all of these samples, 
and the optical parameters of LS are implemented based on 
precise measurements [37-43]. Comprehensive optical pro- 
cesses are simulated using Geant4 [44]. Furthermore, the of- 
ficial electronic simulation (including SPMT’s charge smear- 
ing, transit time spread and dark noise, etc, which are refer- 
enced the measurement in [45]) and charge reconstruction 
are also applied in this study. 

In the calibration, when the ®°Ge source is loaded at one 
of the 227 planned locations, the probability of the unfired 
state for the it” SPMT can be estimated as Eq. 9. Where 
Maz gihreshold corresponds to the number of events with q; < 
ge. and Ntota] is the total number of events for the cal- 
ibration sample, and in this study Ntota: = 10000. According 
to Eq. 6, for convenience, we use an effective mean value of 
detected PEs (2°'-8°""°°) for the it SPMT in Eq. 10. Then 
257 i source can be estimated by Eq. 11 using the calibration 

data. Obviously, this calculation requires that the i” SPMT is 
not fired for all 10000 events in the calibration sample (®°Ge), 
otherwise this method is no longer applicable. Due to the 
small visible energy of °°Ge and the fact that even the most 
marginal of the 227 calibration positions is about 2 meters 
away from its neighboring SPMTs, the above extreme sce- 
nario is extremely unlikely. 
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threshold true_source 
Panfired = PUG < qi |u; ) 
265 B Ng; <qthreshold (9) 
Ntotal 
__,,det_source 
Panfirea = Poisson(k; = 0|u3*20"e*) = e77: 
ia (10) 
det_source 
—_— = — ln Panfirea 
"e Noga (1) 
= -In 
Ntotal 


phy_source 


28 It should be noted that p; 
29 energy reconstruction (Eq. 8), while u; includes the 
270 contributions from visible energy, PMT dark count, charge 
smearing and the threshold effect. According to Eq. 6, Eq. 10 


is the value required for 
det_source 


271 


272 and PMT parameters (98°, 
273 bration, we can find the relationship between yu 


t 
ara putetsource | Fisure 5 shows an example for a SPMT, whose 
275 Spe resolution and dark count rate are 30% and 1 kHz, re- 
27 spectively. The readout window is 1000 ns. It appears that 
277 dark counts dominate p.2¢'-°""°* for small pP?¥S°""°*, while 
272 the combined effect of dark count, smearing and threshold re- 
27 mains stable at around 2% for the given setting as P?Y¥So""ee 
280 increases. 
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Fig. 5. The relationship between p; and u; 
SPMTs whose spe resolution and dark count rate are 30% and 1 kHz, 
respectively; and the readout window is set to 1000 ns. The bottom 


panel shows the ratio of pd°tSource/p,det-source 


Considering that the CD has good symmetry and the 
SPMTs with the same relative position with respect to the cal- 
ibration source exhibit similar responses, in order to enhance 
the accuracy, we further group and combine SPMTs with sim- 
ilar @sprr values in the calculation of u? hy-source Th this 
work, @spmr is divided into 1440 groups from 0° to 180°, 
with 0.125° per group. The same approach has been success- 
fully verified and applied in [23]. As a result, for each °°Ge 
source location P(r, 6,6 = 0), SPMTs in the same ĝspmT 
group have similar values of jP?Y-°°""°, the average of which 
is denoted as „PPY -source (7, Ospmr). The calibration map 
can be constructed after ĉ8Ge scans 227 locations on the X-Z 
plane and all ppPhYS°"rce (7 Aspyr) calculated. Considering 
the calibration performance and time consumption in JUNO, 
there are only about 227 calibration points available in the 
current calibration strategy. So it’s necessary to apply inter- 
polation for the remaining positions. Figure 6(a) and Fig. 6(b) 
23 Show examples of the calibration map before and after inter- 
299 polation, respectively. 
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(b) Ogparr in [18°,18.125°) 


Fig. 6. Example of calibration map for one 0spmr group whose 
pmr Values from 18° to 18.125° 


The following is a brief summary of the main steps in con- 
structing the calibration map. 


(1) For ®8Ge loading at a location P (r, 0, œ = 0), calculate 


the effective mean value of the detected PEs Tae laced for 


each SPMTs using Eq. 11; 


(2) Calculate pPhY-seurce 


and threshold effect; 
(3) Calculate phy source (77 Ospmr), which is the average 
value of phy Sern 


i 


by correcting PMT dark counts 


for each Ogpmrr group; 


(4) Repeat the above steps for calibration data at all loca- 
tions; 

(5) Apply interpolation to the remaining positions; 

(6) For a given ®8Ge location 7 (r,0,¢ 0), the 


pPhy_source(G? @spxrp) value corresponding to each SPMTs 
can be obtained; 


The calibration map is generated using the calibration data 
on the X-Z plane (@ = 0) by considering that the detector 
exhibits good symmetry along the ¢ direction. According to 
a detailed study in [23], the ¢ symmetry is reliable for event 
located within 16 m, while the ¢ dependence can not be ig- 
nored for region near the edge. In the future, the ¢ symmetry 
needs to be checked and validated using real data, and it can 
be corrected if necessary. 
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C. Construction of maximum likelihood function 


To reconstruct the visible energy of a cluster-like event 
whose energy-deposit center is known in the detector, a like- 
lihood function can be constructed as follows: 


Nunfired Ntirea 


L Panfired Parea (12) 


1 1 
where N = Nunfirea + Nirea = 25600, Nunfired and Nirea 
correspond to the number of SPMTs with q; < 0.3 PEs and 
qi = 0.3 PEs, respectively. Panfirea and Pgreq are the prob- 
abilities of unfired state and fired state, which can be calcu- 
lated using Eq. 6, respectively. In the calculation, ??Y of 
each SPMT can be estimated by considering the relationship 
in Eq. 8 and invoking the pP}y-5°0urce (7, Aspyr) value from 
the calibration map. In addition, differences in quantum ef- 
ficiency (QE) between individual SPMTs in the same Ospyrr 
group should be considered. As a result, ?Y-*°""“° needs a 
correction and it can be calculated in the following: 


phy_source __ QE; 


‘ 4 ya QE; 


x pPhy_source (T, Ospmrr) 


(13) 


where Q E; is the QE of the it? SPMT, m is the number of 
SPMTs in the Ospmr group which is classified by the 6spmrr 
value, and QE; is the QE of the jt” SPMT in this group. 
Next, the ROOT’s minimization class TMinuit [46, 47] is 
used for minimizing — ln £. In minimization, the visible 
energy Eyis of the cluster-like event (whose energy-deposit 
center is known and will be introduced in Sec. IV A) is the 
parameter to be determined, and its initial value (initia) 
can be estimated as Eq. 14 using the total number of PEs 
(totalPE) of all SPMTs to reduce the reconstruction time. 
es ) is total PE observed by 25600 SPMTs for 
he ©8Ge source located at different positions (Fig. ??). Com- 
paring the value of total P E5°™°° (r) at the center of CD and 
at around r = 15 m, there is about 7% non-uniformity intro- 
duced by the reception of PMTs and the optical attenuation 
effect in the process of photon transmission. As for the de- 
creasing for radius larger than ~15.6 m, it is mainly caused 
by the total reflection effect and the shadowing effect. 


totalPE(r) 
total P Esouree (r) 


After minimization, TMinuit returns the predicted visible 
energy value as the result of energy reconstruction. Further- 
more, the response difference caused by the spatial scale of 
energy deposition is investigated in Sec. III D and it is found 
that has little influence to our analysis. 


Finitial = Esome’ (14) 


D. Comparison of cluster-like event and point-like event 


68Ge is a positron source. The positron annihilated in 
the source capsule (stainless steel + PTFE) with a pair of 
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404 


0.511 MeV ys emitted. Most of the energy of ys is deposited 
within ~30 cm in the LS, so ©*Ge is not strictly a point-like 
source, but a cluster-like source similar to the high-energy 
electrons described earlier in Sec. I. In this paper, in order 
to determine the location-dependent p??¥°°""* from calibra- 
tion data and carry out energy reconstruction, we construct a 
calibration map using the approximate point-like Ge source 
(cluster size in ~30 cm). In the energy reconstruction of 
high-energy electrons with larger cluster sizes (several me- 
ters), considering that their energies are deposited in a rel- 
atively symmetrical LS space without (or with little) distor- 
tion caused by energy leakage, and the centers of their energy 
depositions similarly coincide with their locations, the cal- 
ibration map is expected to be applicable by assuming that 
all energy deposition is equivalent to occuring at the energy- 
deposit center. However, this assumption and simplified treat- 
ment may bring bias to energy reconstruction, so we need to 
investigate the response difference caused by the spatial scale 
of energy deposition for high-energy event. 

According to Eq. 8, for an event with a known position in 
the LS, with E£S°""°* and aes Source obtained from the cal- 
ibration map, as shown in Eq. 15, the sum of all SPMT’s 
pP” is directly proportional to the visible energy of the event. 
Therefore, any bias presents in the sum of all SPMT’s phy 
indicates the deviation of the visible energy. 


(15) 


3 iH phy _ > iH phy source 


In energy reconstruction, our algorithm assumes that all en- 
ergy deposition is equivalent to occur at the energy-deposit 
center. However, in real detection, ae of the it” SPMTs 
is contributed by the cumulative effect of each secondary en- 
ergy deposition. To estimate the potential bias caused by the 
spatial scale of energy deposition for high-energy event, the 
sum of all SPMT’s u? phy (Eq. 15) is adopted using different 
strategies. As shown in Eq. 16 and Eq. 17, they are calculated 
as point-like event and cluster-like event, respectively. 


Fesource 


N 

> j phy = 
Hi |point-iike = 
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(16) 
N Nag 
h; 
X x` 5 AEw x pe y_ -source (>>) 
i=1 a=1 
N 
5 yr | cluster—like = 
i=1 
N Nae (17) 
Fysource >» 5 AEa x mee source (72) 
t=1 a=1 
N h 
Ratio = Pe MA Y leiuster jike (18) 
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De ie | point—like 


45 The official simulation software of JUNO is applied to gen- 
4% erate electron samples with different kinetic energies, and the 
47 details of energy deposition (AE, and T) are recorded for 
the calculation of Eq. 16 and Eq. 17. Then we compare the 
calculation results as shown in Eq. 18 and Fig. 7. It was found 
that the ratio is close to 1 with the increase of R°, and it has 
a small bias (< 0.5%) at the edge for electrons whose kinetic 
energies are larger than 200 MeV. This result indicates that 
our algorithm and reconstruction strategy is applicable to the 
energy reconstruction of high-energy events which have large 
spatial scales of energy deposition. 
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Fig. 7. Ratio in different positions 


IV. RECONSTRUCTION RESULT 


416 


In this section, the JUNO offline software is used to simu- 
ais late electrons uniformly distributed in the CD with different 
ə kinetic energies (10, 20, 50, 100, 200, 350, 500, 700, 1000, 
and 2000 MeV) as the MC sample to validate the reconstruc- 
tion algorithm. The MC samples are generated after a full- 
chain simulation in JUNO by including realistic detector ge- 
ometry, comprehensive physical processes and optical pro- 
cesses, official electronic simulation and charge reconstruc- 
tion. A 16 m radius cut is applied to avoid energy leakage for 
high-energy event and the total reflection effect at the edge 
of detector. Finally, there are about 10000 events for each 
electron sample. 
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429 A. Reconstruction of energy-deposit center 


As introduced in Sec. IIIC, event’s energy-deposit cen- 
ter is required to reconstruct the visible energy. So, be- 
fore energy reconstruction, it’s necessary to reconstruct the 
energy-deposit center. Compared to the vertex reconstruc- 
434 tion of point-like events, the energy-deposit center recon- 
ass struction of high-energy cluster-like events faces greater dis- 
ass persion (Sec. I). After investigation, it is found that the time- 
437 based algorithm developed and verified by [48] is suitable and 
438 can be applied in our analysis. 


430 


43 


432 
433 


439 The time-based algorithm uses the distribution of time-of- 
o flight (TOF) corrected time At (Eq. 19) of an event to re- 
construct its vertex and to (event time). The principle of the 
time-based algorithm is that the At distribution is indepen- 
dent from the event vertex after applying the time-of-flight 
correction. In this paper, we apply it to the reconstruction of 


the energy-deposit center. 
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446 At; = t; — tofi (19) 

In Eq. 19, t; is the first hit time of the it? SPMT and tof; is 
the time of flight from energy-deposit center to the it” SPMT. 
In the calculation of to f;, the optical path length includes both 
the length in the LS and in the water. Then a correction vector 
will be constructed and minimized by iterating the energy- 
deposit center. More details can be found in [48]. Finally, 
4s3 the reconstructed energy-deposit center can be obtained and 
454 Fig.8 and Fig. 9 show the performance. The reconstruction 
ass biases of 6 (Fig. 8(b)) and @ (Fig. 8(c)) remain small and sta- 
4s6 ble as the energy increases, while the reconstruction bias of R 
(Fig. 8(a)) gradually increases with energy. However, it can 
still be controlled within 150 mm at 1 GeV. Considering that 
the cluster size could be several meters for a 1 GeV electron, 
this bias is still acceptable. And the effect of the deviation 
on energy reconstruction will be discussed later. On the other 
hand, the reconstruction bias of R tends to decrease at the 
edge of the detector compared to other regions, mainly due 
to energy leakage near the edge, especially for high-energy 
events. 

In Fig. 9, it can be found that the reconstruction resolutions 
of R, 0 and ¢ increase with electron energy between 50 MeV 
and 1 GeV. For example, the resolutions of R, 0 and ¢ are 
about 100 mm, 0.5° and 1.0° for 50 MeV electrons, respec- 
tively, while the resolutions of R, 0 and ¢ are about 340 mm, 
1.8° and 3.0° for 1 GeV electrons, respectively. This effect 
is mainly due to the greater dispersion of the energy depo- 
sition of high-energy electrons. The resolutions of 10 MeV 
electrons are slightly larger than the resolutions of 50 MeV 
electron. This is because the hit number of SPMTs is small 
(~400 PEs for 10 MeV) and less information is available for 
the reconstruction using the time-based algorithm. 


447 
448 
449 
450 


45 


452 


457 
458 
459 
460 


46 


462 
463 
464 
465 
466 
467 
468 
469 
470 


47 


472 
473 
474 
475 
476 


477 


478 B. Energy reconstruction performance 


Next, the performance of energy reconstruction will be in- 
troduced. The reconstructed energy spectra for electrons with 
different kinetic energies are shown in Fig. 10. The blue spec- 
tra correspond to events whose energies are fully contained 
(FC) in the LS, while the green spectra correspond to events 
whose energies are partially contained (PC) in the LS. For 
electrons with energies greater than 500 MeV, the proportion 
of PC events becomes larger and the 16 m cut can not to- 
tally excluded the case of energy leakage. The FC spectra 
can be well-fitted with a Gaussian function and the recon- 
structed energy is about 6% larger than the deposited energy 
of electron. According to the official simulation in JUNO, 
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Fig. 9. Reconstruction resolution of energy-deposit center 
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when anchored at the 2.223 MeV gamma peak generated by s2 reconstruction are applied in the reconstruction result shown 
42 (n, Y)H, the high-energy electron has an energy non-linearity 52: by the pink stars and the green triangles. Comparing energy 
43 Of ~6% [49]. Thus, this deviation is understood and it is s2 resolutions in different conditions, it can be found that the 
44 Mainly caused by the energy non-linearity response of LS. s23 energy reconstruction performance of the high-energy events 
49s On the other hand, the non-uniformity of energy reconstruc- 524 is good using the OCCUPANCY strategy and the energy res- 
4s tion may also introduce some small deviations, but generally s2 olution is about 0.8% for 1 GeV electrons in the ideal case 
497 less than 1%, which is shown in Fig.11(c). see (red solid squares). For a more realistic situation by includ- 
48 To understand the non-uniformity shown in Fig.11(c), 527 ing electronic simulation and charge reconstruction, the res- 
499 Fig.l 1(a) and Fig.11(b) can be compared, which corre- sz olution is only about 0.3% worse which indicates that the 
so sponds to the cases using true energy-deposit center with- s2 correction works well in controlling the influence from PMT 
so out/with electronic simulation and charge reconstruction. In ss dark count and threshold effect. From the comparison of 
so Fig.11(a), for electron samples with different energies, the ss: solid points and hollow points, the PC events mainly affect 
so non-uniformity is consistent at about 0.5% from the center of s32 the electrons whose kinetic energy larger than 100 MeV and 
so detector to the edge. After electronics simulation and charge sss their energy resolutions will deteriorate by about 1%. In real 
reconstruction (Fig.11(b)), there is a slight increase in non- s34 detection, the reconstructed energy-deposit center is required 
uniformity, but it still remains within 1.5%. Fig.11(c) cor- sss for energy reconstruction, and its smearing will introduce ad- 
responds to the case using reconstructed energy-deposit cen- s3 ditional smearing on the reconstructed energy, especially for 
ter which shows deviation (Fig.8), as a result, for 500 MeV 37 the high-energy electrons. As a result, the energy resolution 
so and | GeV electrons, their energy non-uniform are about 2% sss is about 3.2% for 1 GeV electrons based on our algorithm. In 
sto and 3% at the edge. Furthermore, if the PC events are in- s3 Fig. 14, the relationship of energy resolution and the fired ra- 
cluded (Fig. 12), they mainly affect the non-uniformity of s« tio of SPMT is investigated using the electron sample without 
st2 high-energy electrons located in the edge region. s41 electronics simulation, and true energy-deposit center is ap- 
Figure 13 shows the performance of energy resolution. The s42 plied. In general, a higher fired ratio of SPMT corresponds to 
514 SOlid points correspond to the reconstruction results of FC s« a better energy resolution. This indicates that our algorithm 
events, while the hollow points include both FC events and s4 has the potential to be applied to higher energy events when 
PC events. The red squares and pink stars denote the cases sas the fired ratio of SPMT is not close to 1, but must solve the 
517 using the true energy-deposit center for energy reconstruc- s4 problem of energy-deposit center reconstruction which has a 
tion, while the green triangles use the reconstructed energy- s47 larger bias at higher energies. 

st9 deposit center. In addition, electronics simulation and charge 
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Fig. 10. Discrete energy reconstruction with reconstructed edep vertex after electronic simulation and charge reconstruction. The blue line is 
the FC events and green line is the PC events. According to the fitting results (red line) of FC spectra, it can be observed that the reconstructed 
visible energy is about 6% larger than the deposited energy of electron. More specifically, for electrons with kinetic energies of 10 MeV, 
50 MeV, 200 MeV, 500 MeV, 700 MeV and 1 GeV, the ratio of reconstructed visible energy to deposited energy is found to be 1.062, 1.061, 
1.058, 1.059, 1.060 and 1.064, respectively. The corresponding explanation is provided in the text 
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Fig. 11. Uniformity of discrete energy reconstruction for FC events. On each plot, black vertical dotted lines corresponds to R = 15 m, 
and red vertical dotted lines corresponds to R = 15.6 m, which is the boundary of the total reflection region (caused by the refractive index 
mis-match between LS and water) 


Accurate energy reconstruction is crucial to detect various 
physics events in a wide energy range from MeV to GeV in 
JUNO. This work focuses on the energy reconstruction of 
sub-GeV event based on 3-inch PMTs and the OCCUPANCY 
method. Our reconstruction shows good performances in the 
test of MC simulation samples. The energy non-uniformity 
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can be controlled within 1% from the center of detector to 
the edge for electrons whose kinetic energies are smaller than 
500 MeV. As for 1 GeV electron, the energy non-uniformity 
can be controlled within 3%. 
1 GeV electron FC events and FC+PC events are about 2.7% 
and 3.2%, respectively. Our algorithm has the advantage of 
small dependency on precise charge reconstruction by mainly 
using the information from SPMT’s firing states. 
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(c) With electronic simulation and charge 
reconstruction, using reconstructed energy-deposit 
center for energy reconstruction 


On each plot, black vertical dotted lines corresponds to 


R = 15 m, and red vertical dotted lines corresponds to R = 15.6 m, which is the boundary of the total reflection region (caused by the 
refractive index mis-match between LS and water) 
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Fig. 13. Resolution in several discrete energy in simulation phase. 
The red squares correspond to the resolution in detector simulation. 
The pink stars and the green triangles correspond to the resolution 
after electronic simulation and charge reconstruction. The energy re- 
constructions of the red squares and pink stars use the true center of 
energy deposition, while the green triangles utilize the reconstructed 
center for their energy reconstruction. The energy reconstruction of 
the red squares and the pink stars using the true energy-deposit cen- 
ter, while the energy reconstruction of the green triangles uses the 


reconstructed energy-deposit center. 
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